TTP 02-13 
UPRF-2002-09 



Two-loop amplitudes with nested sums: Fermionic 
contributions to e^e~ qqg 



Sven Moch'^, Peter Uwer*^ and Stefan Weinzierl' 



"Institut fiir Theoretische Teilchenphysik, Universitdt Karlsruhe 
76128 Karlsruhe, Germany 

^ Dipartimento di Fisica, Universita di Parma, 
INFN Gruppo Collegato di Parma, 43100 Parma, Italy 



Abstract 



We present the calculation of the n y -contributions to the two-loop amplitude for e^e^ qqg 
and give results for the full one-loop amplitude to order £^ in the dimensional regularization 
parameter. Our results agree with those recently obtained by Gai^land et al.. The calculation 
makes extensive use of an efficient method based on nested sums to calculate two-loop inte- 
grals with arbitrary powers of the propagators. The use of nested sums leads in a natural way 
to multiple poly logarithms with simple arguments, which allow a sti'aightforwai^d analytic 
continuation. 



1 Introduction 



The last thirty years of experimental studies at colliders together with the theoretical investi- 
gations, have taught us that perturbative quantum chromodynamics (QCD) gives an excellent 
description of short-distance scattering of strongly interacting partons. Indeed, today the theory 
has reached a sufficient maturity that it is no longer the target of experimental studies, but rather 
a tool in the search for new physics beyond the standard model. 

The search for new physics in particle physics, being pursued at present and upcoming collider 
experiments at the Tevatron and the LHC, rely on our ability to make precise predictions for 
QCD and QCD-associated processes. The accuracy reached already in present collider experi- 
ments demands next-to-next-to-leading (NNLO) theoretical predictions within the framework of 
perturbation theory. For example the strong coupling constant a^, whose precise value affects 
many cross sections, can be measured by using the data for e+e^ — > 3 jets. At present, the er- 
ror on the extraction of a, from this measurement is dominated by theoretical uncertainties [jT]], 
among the main sources there being the truncation of the perturbative expansion at a fixed or- 
der. Up to now, event shapes in 3-jet events have been calculated at next-to-leading order (NLO) 
for massless and massive quarks [f7|]-[[TTl]. To reduce the theoretical uncertainties, it is 

necesarry to extend the calculation for massless quarks to next-to-next-to-leading order. The 
calculation of e^e^ 3 jets at NNLO requires the tree-level amplitudes for e^e~ — > 5 partons 
inp, the one-loop amplitudes for e+e^ 4 partons [P^-[P^ as well as the two-loop am- 
plitude for e^e~ qqg together with the one-loop amplitude e+e qqg to order £^ in the 
dimensional regularization parameter. 

While for inclusive quantities like the total hadronic cross section in e+e" -annihilation even 



higher orders have been calculated in the past [^SQ-PQ], the calculation of two-loop four-point 



scattering amplitudes has been the main obstacle for a long time. Due to tremendous activity in 
that field during the past three years this problem can be considered to be solved - at least 
for the case of massless internal quarks and only one external massive leg. 

The by now more or less standard approach to calculate two-loop four-point functions has been 
inspired by the techniques developed in the calculation of two-point functions. The starting point 



is a reduction of the tensor integrals through Schwinger parametrization [P2|]-[p4|]. This yields 
immediately scalar integrals with a higher dimension and raised powers of the propagators. The 
usual way to proceed then, is to apply repeatedly algebraic relations between these integrals, 
which follow from Poincare- and Lorentz-invariance [^-[^. Finally one ends up with a small 
set of so-called master integrals which must be solved analytically. While for simple topologies 
it is often straightforward to find the reduction scheme (i.e. 'triangle rule'), in general it is a 
non-trivial task to solve this problem. 



In a recent publication [ P8| , p9[ |, we have therefore proposed a different method to attack this 
problem. The basic idea is the following. We solve the scalar integrals in higher dimensions and 
with raised powers of propagators directly in terms of nested sums instead of reducing all the 
integrals to a small set of master integrals. The aim of this paper is to illustrate this method in 
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the calculation of the fermionic contributions to the two-loop amplitude e+e qqg, i.e. the 
contributions proportional to the number of quark flavours nf. We present our results in terms 
of multiple polylogarithms [pO|]-[^, which arise naturally from the use of nested sums. In 
addition, we show that these multiple polylogarithms can easily be continued analytically. As 
a consequence the amplitudes for e+e^ — > qqg presented here can also be used for (2 + l)-jet 
production in deep-inelastic scattering and the production of a vector boson (V7, Z or Drell-Yan 
pair) in hadron-hadron collisions. The respective amplitudes can be obtained by the crossing 
symmetry and simple coupling constant modifications. 

The outline of the paper is as follows. In section ^ we present a few properties of the two-loop 
amplitude. In particular, we discuss the kinematics, the ultraviolet (UV) renormalization, and 
the structure of the soft and coUinear singularities. In section ^ we outline the calculation. In 
the following section ^ we present our results and compare them with those recently obtained by 
Garland et al. . We give our conclusions in section ^. Appendix contains the results for the 
one-loop amplitude to order £^ and appendix^ summarizes properties of multiple polylogarithms 
under analytic continuation. 



2 Preliminaries 



2.1 Kinematics 



In the following we study the reaction 

e+ ^ <? + g + ^. (1) 
We treat all quarks as massless, that means we work in QCD with nf massless quark flavours. 



To be consistent with earlier work 014] , |15| ] we calculate the amplitude for the reaction with all 
particles in the final state 



0^q{pi)+g{p2)+q{p3)+e {p4)+e+{p5). (2) 
The kinematical invariants are denoted by 

Sij = [pi + PjY, Sijk = [pi + Pj + Pkf , 5 = 5123 , (3) 

and it is convenient to introduce the dimensionless quantities 

Xi = , X2 = . (4) 

■^123 ■^123 

For pure photon exchange, the complete amplitude Jij for e+e^ qqg can be written as the 
product of a lep tonic current with the hadronic current H^: 

= -'-e^Q''Ti;L,H^' = -e^Q'^nAy, (5) 
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with denoting the electric charge of the outgoing quarks in units of the elementary charge 
e = -\/47roc. The generator of the SU(A^) gauge group is given by T'^. The indices ij and a 
describe the color of the outgoing quarks and gluon. The normalization of the color matrices is 
taken to be 

Tr(^r"r*) = l5°^ (6) 

The leptonic current is given by 

Lfj = u{p4)y^iV{p5), (7) 

with u, V denoting the spinors of the outgoing leptons. As we will show later, it is sufficient to 
consider pure photon exchange: Working in a helicity basis, the pure photon exchange amplitude 
j^y allows the reconstruction of the full amplitude including Z-boson exchange by adjusting the 
couplings. Using the anti-commutation relations for the y-matrices, one can always achieve the 
following form of the hadronic current: 

H^, = Ci^{pi\^2\P3)£g^ + C2^{pi\h\P3){£g-Pl)P3^j 

+ C3^{pi \h\P3) (Eg • P3)Pl^, + C4^{pi\^2\P3){£g ■ P\)Plt, 

+ C5 {pi \h\P3) (Eg ■ P3)P3f, + (pi IY^Ips) (Eg ■ p\) 

+ Ci-^{P\\l^P3){^g- P3) +Ci-^{P\\4g\P3)Plfji 

+ C9\p\\4g\P3)P3t, + CiQ\pi\4g^2y^,\P3) 

+ Cii ^ {pi \ig\p3)P2^, + Cn-^{p\\h\P3){^ ■ Pl)P2fj 

+ Ci3-^{pi\h\P3)i£g- P3)P2^, (8) 

where we have used {pi \ respective \p3) as a short-hand notation for the spinors u{pi) and v(p3) 
of the outgoing quark and anti-quark. The dimensionless functions c, depend only on the ratios 
Xi, the spacetime dimension d = 4 — 2e and the renormalization scale /j, 

Ci = Ci(xi,X2,£,^^ . (9) 
Due to various constraints, for example current conservation, 

iPl, + P2, + P3,)H^ = 0, (10) 

the functions are not all independent of each other. It can be easily chown that C2,C4,C6,ci2 
are sufficient to reconstruct all remaining functions. A similar conclusion has been drawn in ref. 
[p3|]. The relations between C2, C4, ce, cn and the remaing functions are: 

C3(xi,X2) = -C2{x2,Xi), 
C5{xuX2) = -C4{x2,Xi), 
Cl3(-^l,-^2) = -Cl2(^2,^l), (11) 
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and 

1/1 N 1/1 N ^1 1 / \ 

Cl = -:z{^-X2)C3- -{l-Xi)C5 + —C6- -{xi+X2)Ci3, 
2 2 X2 2 

Xl 

Cl = C6, 

X2 

1 1 

1 1 

eg = -2-^iC2-2-^2C5, 

1/1 1/1 1 (-^1 +■'^2) 

ClO = +-{l-Xi){C2-C5)--{l-X2){C3-C4)+- Ce 

4 4 2 X2 

+ +X2)(ci2-Ci3), 

Cll = l;{l-X2)c3 + ]-{\-Xi)c5-—C6 + lrXiCi3-l-XiCi2. (12) 
2 2 JC2 2 2 

In the actual calculation we have not used these constraints. Instead we calculated all ci-cn and 
used eqs. ( [TT| ) and (|12|) as a cross-check on our calculation. 

Beyond the leading -order, one encounters UV as well as soft and coUinear singularities. We 



use dimensional regularization [g5], Q to regulate both types of singularities. There are several 
variants of dimensional regularization which are used in loop calculations in QCD: Conven- 
tional dimensional regularization (CDR) [P5|] continues all momenta and all polarization vectors 



to d dimensions. The 't Hooft-Veltman (HV) scheme [ |250 takes the momenta and the helici- 
ties of the unobserved particles in d dimensions, whereas the momenta and the helicities of the 
observed particles are four-dimensional. The CDR scheme is often employed within the interfer- 
ence method, but is not suited for the calculation of amplitudes. Enforcing the CDR scheme in 
the calculation of amplitudes requires the introduction of external states with "e"-helicities [|36|]. 
Furthermore, there are several versions of four-dimensional schemes on the market. Despite the 
name "four-dimensional schemes", they are variants of dimensional regularization. The name 
refers to how these schemes treat unobserved internal particles and the Dirac algebra. The four- 
dimensional helicity scheme (FDH) [BTp introduces an additional parameter ds for unobserved 



internal states, which is set to 4 at the end of the calculation. It has the advantage that it respects 
supersymmetric Ward identities up to two loops. The four-dimensional scheme defined in [pBp 



keeps the Dirac algebra in four dimensions and allows the use of four-dimensional Fierz- and 
Schouten identities, at the expense of having to restore Ward identities. These schemes can lead 
to considerable simplifications, in particular if many external particles are involved. For the pro- 
cess e+e^ — > 3 jets, the number of external particles is relatively small and we do not consider 
these schemes further. 

In this paper we keep our calculation rather general and we decompose the amplitude into spinor 
strings, which we can compute without any reference to the dimensionality of the external polar- 
ization vectors. Internal particles and the Dirac algebra are treated in d dimensions. From these 
spinor strings we can easily deduce the results in the CDR scheme and the HV scheme. The 
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result in the CDR scheme is obtained by interfering the two-loop amplitude with the Bom am- 
plitude. On the other hand, by contracting with explicit representations of polarization vectors, 
we obtain helicity amplitudes in the 't Hooft-Veltman scheme. 

Including the Z-boson exchange, the situation becomes slightly more complicated. Due to the 
presence of 75, a specific scheme has to be chosen. Schemes which allow for a consistent treat- 



ment of Ys are for example the HV scheme or the scheme defined in [pBO. In both schemes, 
the regularization procedure violates certain Ward identities, which have to be restored by finite 
renormalizations. However, since the amplitudes considered in this paper do not contain closed 
fermion loops with axial-vector couplings, the results for Z-boson exchange can be obtained 
from the ones for pure photon exchange by a simple adjustment of the electro-weak couplings. 

It should be noted that in general the finite part of the amplitudes are also scheme-dependent. As 
long as soft and coUinear singularities are considered, this is not really an issue. Using the same 
scheme in the calculation of the divergent contributions from the real corrections any scheme 
dependence cancels out at the end. For the UV divergences however, one has to keep in mind 
that in general the coupling constant is scheme dependent. The values of the coupling constant 
in two different schemes are related by a finite renormalization. The coupling constants in the 
CDR scheme and the HV scheme are identical, since the internal states in these two schemes are 
treated in the same way. The CDR or HV scheme, together with the modified minimal subtraction 
prescription (MS), defines the usual coupling and it is most useful to quote our results in 
these schemes. 

To obtain the results in the 't Hooft-Veltman scheme we work in the helicity basis. For fermions, 
spinors of definite helicity are given by: 

"±(p) = ^(i±Y5) 

v±ip) = ^{ItJsHp). (13) 
We introduce the following short-hand notation for spinors with definite helicity: 

|z±) = \pi±) = U±{pi) = V^{pi), 

{i±\ = {pi±\ = u±{pi) = v^ipi), (14) 
and the spinor products are then defined as 

{pq) = {p-\q+), 

[pq] = {p+\q~). (15) 
Gluon polarization vectors are expressed in terms of Weyl spinors as [p9|]-[p4|] 
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Figure 1 : Labelling of the external legs for the amplitude. The lepton pair is denoted by p4 and 

P5- 



where k is the gluon momentum and q is an arbitrary null reference momentum. The dependence 
on the reference momentum drops out in final gauge-invariant amplitudes. An appropriate choice 
of q can lead to a significant reduction in the number of diagrams which need to be evaluated. 

Helicity conservation for massless quarks and leptons ensures that there are only 2-^ = 8 possible 
helicity configurations for e^e^ qqg, namely two choices for each fermion line together with 
two possible gluon polarizations. Because the electron line couples through the current (4 ± 
|y |5±) = (5 =F |y |4=f)5 it is trivial to reverse its helicity simply by exchanging p4 ^ and 
by adjusting the weak couplings. Parity and charge conjugation can be used to further reduce 
the number of helicity amplitudes which need to be calculated. Parity reverses all helicities 
simultaneously and is implemented by complex conjugating all spinor products (e.g. {ij) ^ [ji])- 
Charge conjugation reverses the arrows of each fermion line. In addition there is a factor (—1) 
for each external gauge boson. Thus we are left with just one independent helicity amplitude, 
which we take tobeAy(l + ,2+,3",4+,5"). 

Keeping the reference momentum q in the gluon polarization vector arbitrary we obtain: 

A,(l+,2+,3-,4+,5-) = -^^-i{2ci(^5)[42][12](23)+C2^(^l)[12][43](35)[12](23) 

+C3^(^3)[32][41](15)[12](23)+C4^(^1)[12][41](15)[12](23) 

+C5^ (^3) [32] [43] (35) [12] (23) + 2c6(^l) [12] [41] (35) 

+2c7 (^3) [32] [41] (35) + 2c8(<?3) [12] [41] (15) 
+2c9(^3) [12] [43] (35) + 4cio(^2) [12] [42] (35) 

+2cii (^3) [12] [42] (25) + cn- (^1) [12] [42] (25) [12] (23) 

+ci3^(^3)[32][42](25)[12](23)|. (17) 

Using q = P3 and the constraints of eq. (|12|) this can be simplified to: 
Ay(1 + ,2+,3-,4+,5-) = 
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[12] 



-(31)[12] 



.(35) [42] 



(1 -Xi) ( C2 + — C6-C12 J +(1 -X2) (C4-C12)+2C12 



[43] (35) ( C2 + -C6 - C12 ) + [41] (15) (c4 - C12) 



(18) 



We see, that for this choice of helicities and q the amplitude Ay is described by three linear 
combinations of C2, C4, C6, C12. 

The perturbative expansion of the functions c; and Ay is defined through 

c, = V^(.iO..(|),a.,(|)%..,o,„3)), 

A, = V^(Af+(^)A<"+(^)V + 0(a.')). (W) 



Inserting the leading-order result for the c; 



^(0) _ ^(0) _ ^(0) _ ^ (0) 



one gets the following result for the tree amplitude 

aJ°^(1+,2+,3-,4+,5-) = 2V2/- 



(35)^ 



(20) 



(21) 



^ ' (12)(23)(45)' 

In the helicity basis one can easily account for the Z-boson exchange by adjusting the couplings: 

j?yz(l''S2''2,3^3^4^,5''5) = e2(_Q9 + v^^v^^2,2(^))^^(l'^i^2^2,3''3,4^,5''5)^ (22) 

with 



'Ez{s) 



s — ni^ + inizTz 



(23) 



Here, mz and Tz are the mass and the width of the Z-boson. The left- and right handed couplings 
of fermions to the Z-boson are 



f l{-Qfsm^%w f -QfsinQw 

^- = • Q — ^+ = — —a — ' 

sin Qw cos Qw cos Qw 



(24) 



f 

where Qf and are the charge and the third component of the weak isospin of the fermion / 
and Qw is the Weinberg angle. 



2.2 Ultraviolet renormalization 

The amplitudes we present are the renormalized ones, i.e. the ultraviolet subtraction has been 
performed. To obtain the renormalized amplitudes in the MS scheme, one replaces the bare 
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coupling ao with the renormalized coupling (Xs{^i^) evaluated at the renormalization scale /u^ 



ao 



where 



{Anfe 



(25) 



(26) 



is the typical phase-space volume factor in J = 4 — 2£ dimensions, Je is Euler's constant, and (3o 
and Pi are the first two coeffcients of the QCD [3-function: 



with the color factors 
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2 „ 17 2 5 

Ca - :^TRnf, pi = — Q - -CATRHf - CfTRHf, 



17 
~6 



Ca=N, Cf 



2N 



Tr 



(27) 



(28) 



It is convenient here and for the subsequent discussion of the soft and coUinear singularities to 
introduce a different notation In an orthogonal basis of unit vectors \i,i,a) in the three 

parton color space we define an abstract vector | M ) through 



Jl = {ij,a\9\{), 
with the expansion in the coupling a, defined by 



(29) 



(30) 



Then, the renormalized two-loop amplitude can be expessed as 



Kren) 



47ra\/47ra,5, 



-1/2 



At two loops the relation between the renormalized and the bare amplitudes is given by 



ren / 


-"^'bare / ' 


ren / 


= s,^ 


Wale 


ren / 


= 


^b£e 



-1 



i)\ , ^3p2 



bare ^ 



+ 



8e2 4e 



Thus, we obtain for the renormalized functions c,>en 

(l),ren 



c-l„(l)>bare 
-i — ^£ 5 

J2),ren _ 2^(2), bare _ 3p() ^- l (l), bare 

■i ~ 2e ^ ' ' 



(0)\ 
bare / ' 



(31) 



(32) 



(33) 



for / = {2,4, 12} and 



(l),ren _ 1 (l),bare _ Po (0) 

^6 ~ ' 



(2),ren ^ 2 (2),bare _ 3pO^-l (l),bare , _ (2),bm-e ^4) 

6 £6 2z ^ ^ V 8e2 4£ ' ^ ■ 



In this paper we set the renormalization scale iJ^ = s. The complete scale dependence is easily 
recovered by expanding the prefactor 



2\ 



s 

accompanying an Z-loop amplitude to the appropriate order in £. 



(35) 



2.3 Infrared structure 



Based on universal properties of soft and coUinear limits, the infrared pole structure of two-loop 



amplitudes has been predicted by Catani [p3|]. Here, we briefly review how to organize these 
infrared poles for e+e qqg. We start with the one-loop amplitude, which can be written as 



9W' 



l(l)(£) 



(1) 



(36) 



Here l'^^)(£) contains all infrared double and single poles in l/£ and | J '^^)) is a finite remainder. 
At two-loops, the corresponding formula reads: 



(2) 



l(l)(£) 5lf(l)\+l(2)(£ 



The one-loop insertion operator I*^^) is given by 



(2) 



2 \ e 



(37) 



(38) 



where 



and the coefficients and y, are 



1 1 

£^ £ 



T2 rwyZ rw^l /"< 

q — ^q— ^F, ig — t-A, 

3 

yq = yq = :^Cf, yg = ^o. 



(39) 



(40) 
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In general, the color operators TjTy give rise to color correlations. However, the color structure 
for the amplitude e+e qqg is rather trivial and the color operators are proportional to the 
identity matrix in color space: 



T,T, = T,T^ = - 

Explicitly, the one-loop insertion operator reads for e+e 

.2 - £ 



-TrN. 



1 1 



2r(l-£) 



X 



^1 T, 



1 3 1 



2e 



Po 3\ 1 



Ca 2 / £ 



The two-loop insertion operator has the form 

l(2)(£) = 

-il(i)(£) (^l(i)(£) + 2|3o^)+."^^^ 



r(i-2£) 1 



where 



K 



67 



10. 



(41) 



(42) 



Pol+i^^l(i)(2e)+H(2), (43) 



(44) 



The function H^^) is process- and scheme-dependent and for e+e qqg, it is given by p7| ] 



H(2) 



1 1 



1 



.e^,._,H(2) (2) (2) 

4r(i-£)^ £^**'^ ^^^'^ 



(45) 



where H, 



(2) 



H?^ and 



„(2) (T^ 11 2 "tOSA.,! 1^ "ll A 3^ 3 1 

= i4^^-%''+8sj'^+i-4^'-%-108j + i-2^' + T-32j 
25 \ _ 1 



^ Us 2167 ^ 

f 1 fin 

Using the above results we define the finite functions c\ 



l(i)(£)4^^'"" (47) 
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for / = {2,4, 12} and 

(l),fin _ (l),ren T(l)^pN (0) 

^(2),fin ^ 4i)^-_i(i)(e)4i)'-_i(2)(e),(0). (48) 

f2l fin fll fin 

Explicit results for the functions c- " are given in section ^ and for c) in the appendix ^ 



3 Method of calculation 



In this section, we will discuss the method to calculate the virtual amplitudes for e+e qqg. 
We have used QGRAF [ |48| ] for the generation of all Feynman diagrams, which contribute to the 
process e+e^ qqg up to two loops. 

The evaluation of the diagrams leads to tensor integrals, which multiply the various spinor struc- 
tures of eq. Note, that there is no need here to consider projectors for the various spinor 
coefficients. 

The tensor integrals are mapped to combinations of scalar integrals with higher powers of propa- 



gators and different values of d [p2|]-[p4|| . For this purpose, one introduces Schwinger parameters 
for the propagators 



1 1 



:-By r(v) 



dxx'^Uxpixk^). (49) 



Combining the exponentials arising from different propagators one obtains a quadratic form in 
the loop momenta. For instance, for a given two-loop integral with loop momenta ki and k2, one 
has then 



/(J,Vi,...,V^) 




d''k2 1 



m^/2(-yt2)Vl...(_yt2)v„ 

d''k2 



1^ Of^r'^-r' exp E^,*? . (50) 



and 



Y^Xikf = akj + bkl+2cki-k2 + 2d-ki +2e-k2 + f . (51) 

The momenta kj,. . .,kk are linear combinations of the loop momenta ki , k2 and the external 
momenta. The coefficients a, b, c, d^, and / are directly readable from the actual graph: 
a{b) = Y^Xi, where the sum runs over the legs in the ki (^2) loop, and c = Y,^i with the sum 
running over the legs common to both loops. With a suitable change of variables for the loop 
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Figure 2: The PentaBox represents a typical diagram where the triangle relation can be applied. 
In the present calculation, the lower left leg is massive. 



momenta ki , k2, one can diagonalize the quadratic form and the momentum integration can be 
performed as Gaussian integrals over the shifted loop momenta according to 

/^«P(-*^) - (52) 

Lorentz invariance allows immediately to relate the following symmetric tensor integrals to scalar 
integrals: 



and the generalization to arbitrary higher tensor structures is obvious. 

In the remaining Schwinger parameter integrals, the tensor integrals introduce additional factors 
of the parameters Xi and of \/t . These additional factors can be absorbed into scalar integrals 
with higher powers of propagators and shifted dimensions, by introducing operators i+, which 
raise the power of propagator i by one, or an operator d+ that increases the dimension by two. 



1 1 



V,l ' r-T7T = V/ 



— — - j dxi X-' ^ Xi exp {xil<? 



1+ 



At this stage, we are left with sets of scalar integrals of a given topology in 4 + 2m — 2e- 
dimensions, and with raised powers of propagators, where m is a non-negative integer. Then, for 
some topologies like the PentaBox in fig. @, immediate simplifications are possible. By means 
of integration-by-parts [|25|, ^ these topologies reduce to simpler ones. For the PentaBox in fig. 
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Figure 3: The generic CBox. For the CB0X2, the upper right external leg is massive. 



@, for instance, partial integration provides the following triangle relations, 

[(J-2v2-V3-V5)-V33+2"-V55+(2"-l")]PentaBox(m,Vi,...,V7) = 0, 
[(J-V2-2v3-V5)-V22+3"-V55+(3"-4")]PentaBox(m,Vi,...,V7) = 0. (55) 

These relations can be used to eliminate the propagators 1 and 4, such that the PentaBox becomes 
reducible to the CB0X2 shown in fig. (Q). 

After having performed obvious simplifications based on triangle relations, we are then able, to 
calculate directly all necessary scalar integrals with the method of nested sums pE|]. For each 
topology, this requires analytical solutions valid in any dimension and for any (not necessarily 
integer) power of the propagators in terms of nested sums. For the fermionic contributions to 
e+e^ — > qqg up to two loops, it suffices to have these analytical expressions for the CB0X2 pS|], 
the one-loop box with one external mass [ p9[ ] and the one-loop triangle with two external masses 

We give the explicit representations as nested sums for the basic integrals. As a short-hand 
notation, we use V/y = +Vj for sums of powers of propagators in the following. The one-loop 
triangle with two external masses is defined by 

Tri2(m,Vi,V2,V3;xi) = (-.123)"™+'+"'^^ / "^^.^ . Iv, . iv, . Iv. . (56) 



where k2 = ki — pi — p2, k^, = k2 — P3. It can be written as a combination of hypergeometric 
functions 2^1 • The series representation for this integral is given by 

r(£-m + V23)r(l-£ + m-V23)r(m-£-Vi3) 
Tri2 m,Vi,V2,V3;.^l = t-a, ^t-^, ^r^, ^ro ^ (^^^ 

r(vi)r(v2)r(v3)r(2m - 2£ - V123) 



^ z' 



■ nr-£-V23 + - £ -f m - V2) r(/i+V3)r(/i-m + £+Vi23) 



r(/i + 1 + m - £ - V23) r(zi + 1 - m + £ + V23; 

The one-loop box integral is defined by 

Box(m,Vi,V2,V3,V4;jci,JC2) = (58) 

d^'ki 1 



^ ^''^ J {-kjY' {-kiy {-kiy {-kj) 



V4 ' 
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where ^2 = — Pl> ^3 = ^2 — P2 and ^4 = A:3 — ^3. This integral can be expressed in terms of a 
combination of Appell functions of the second kind and the series representation is given by: 



Box(m,Vi,V2,V3,V4;xi,X2) = 

r(/r;-e-Vi23)r(/r;-e-V234)r(l +V123 -/?7 + e)r(l +V234 -/?7- 



r(vi)r(v2)r(v3)r(v4)r(2m - 2e - V1234) 



x; Tiix +V3)r(/2+V2)r(/i + J2 -m + e + Vi234) 



r(/i + l - m + e + Vi23)r(/2 + l - m + e + V234) 



m-e-vi23 



r(ji + m - e - vi2)r(j2 + V2)r(/i + /2 + V4) 



r(/i + 1 + m - £ - V 123 )r(/2 + 1 - m + e + V234) 
^^-e-v234 r(/i + V3)r(i2 + m - e - V34)r(/i + /2 + Vl) 
2 r(ii + l - m + £+Vi23)r(/2 + l + m-£-V234) 



_^^OT_e-vi23^m-e-V234 r(/i + w - £ - Vi2)r(?2 +ot - £ - V34)r(M + ^2 +^ - £ - V23) 
^ ^ r(«i + l+m-£-Vi23)r(i2+l+W-£-V234) 

Finally, the two-loop CB0X2 is defined by 

CBox2(m,Vi,V2,V3,V4,V5;xi,X2) = (60) 

r_,„,r2-+2e+vi2345 f^_^ I I I L_ 

^ J indn J (_;fc2)Vi ^_^2)V2 (_^2)V3 ^_^2)V4 (_^2)V5 ' 

with h = ki + 15 — pi, — I2 — p2, k4 — ki —pwi- The formula for this integral is given by: 

CB0X2(m,Vi,V2,V3,V4,V5;xi,X2) = 

r(2m - 2£ - Vi235)r( 1 + V1235 - 2m + lt)T{lm - 2£ - V2345)r( 1 + V2345 - 2m + Iz) 

r(vi)r(v2)r(v3)r(v4)r(v5)r(3m - 3£ - V12345) 

r(m-£-V5)r(m-£-V23) y y x| 
r(2m-2£-V235) hh'^-'^- 



X 



r(/i + V3)r(/2 + V2)r(/i + i2 -Im + lt + Vi2345)r(/l + /2 - m + £ + V235) 

r( ii + 1 - 2m + 2£ + V1235 )r( j2 + 1 - 2m + 2£ + V2345 )r( ii + ^2 + V23 ) 

2m-2e-Vi235 



r(/i + 2m - 2£ - Vi25)r(/2 + V2)r(/i + /2 + V4)r(/i + /2 + W - e - Vl) 

r(ji + 1 + 2m - 2£ - Vi235)r(/2 + 1 - 2m + 2£ + V2345)r(/i + i2 + 2m-2e- V15) 

2m-2E-V2345 

•*2 

r(/i + V3)r(/2 + 2m - 2£ - V345)r(ii + ;2 + Vi)r(/i + /2 + m - £ - V4) 
r(/i + 1 - 2m + 28 + Vi235)r(/2 + 1 + 2m - 2e - V2345)r(ii + i2 + 2m - 2e - V45) 
2m-2e-vi235 2m-2e-v2345 ^(/i + 2m - 28 - Vi25)r(/2 + 2m - 2e - V345) 

r(/i + l+2m-2e-Vi235)r(/2+l+2m-2e-V2345) 
r(/i + i2 + 2m - 2e - V235)r(/i + (2 + 3m - 3e - V12345) " 
r(/i + /2 + 4m - 4e - V12345 - V5) 



(61) 
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Then, the evaluation of the multiple nested sums proceeds systematically with the help of the 
algorithms of These algorithms rely on the algebraic properties of the so called Z-sums and 
allow to solve by means of recursion the nested sums in terms of a given basis in Z-sums to any 
order in the expansion parameter £. Z-sums are defined by 



Z{n;mi,...,mk;xi,...,Xk) = ^ 

n>ii>;2>...>;A:>0 

They are generalizations of Euler-Zagier sums [BTl p2p 



X, 



-mi,...,mi 



Z{n\mi,...,mk\\,...,\] 



X, 



Ik 



ink 



(62) 



(63) 



or of harmonic sums [|53|] - [p6|] involving multiple ratios of scales. The latter are well known 
from calculations of Mellin moments of deep-inelastic structure functions [|53|], [ |57| ] - [[60|]. An 
important subset of Z-sums are multiple polylogarithms [pO|]-[^, obtained by letting n go to 
infinity in eq. (F^ 



Li 



m^,...,mi 



{xk,...,xi) = Z{oo;mi,...,mk;xi,...,Xk). 



(64) 



The key feature of Z-sums is the fact, that they interpolate between Goncharov's multiple poly- 
logarithms and Euler-Zagier sums, which occur in the expansion of F-functions, 

r(n-F£)=r(i+£)r(n) 

X {l+eZi{n-l)+e^Zn{n-l)+e^Zni{n-l) + ...+e"-^Zn..A{n-l)). (65) 

In summary, the evaluation of the multiple nested sums proceeds by expanding the F-functions 
to the desired order, and by using then the algebraic properties of the Z-sums. In this way, we 
could calculate all loop integrals contributing to the one- and two-loop virtual amplitudes very 
efficiently in terms of multiple polylogarithms. All algorithms for this procedure have been 
implemented on a computer in symbolic manipulation programs. We have chosen to work with 
FORM [§TJ] and in the GiNaC framework §3p. 



4 Results 

Let us now discuss our results for the finite coefficients c-^^'^" through order and the rif- 
contributions to cf'^'^^. 

We have made the following checks on our result. First of all, after UV renormalization, the 



infrared poles agree with the structure predicted by Catani [p3|]. This provides a strong check of 
the complete pole structure of our result. Secondly, having calculated all coefficients C1-C13, we 
could use the various relations eqs. ( [TT] ) and ( [T2| ) between C1-C13 as a cross-check. Note that for 
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C6 there is an additional symmetry: one can show that the combination xiC(, is symmetric under 
the exchange {xi ^ X2), 

xiC(,{xi,X2) = X2Ce{x2,xi). (66) 

Finally, we could compare with recent work of Garland et al.. They have obtained by an indepen- 
dent method the result for the squared matrix elements [B^, i.e. the interference of the two-loop 



amplitude with the Bom amplitude, and the interference of the one-loop amplitude with itself, as 
well as the result for the one- and two-loop amplitude p3|]. Their results are given in terms of 
one- and two-dimensional harmonic polylogarithms, which form a subset of the multiple poly- 
logarithms [pOp-Pl]. Thus, we have performed the comparision analytically. We agree with both 



of their results. 

Here, we prese 
amplitude are deferred to appendix 



f2l fin 

Here, we present explicitly the function c- / = 2,4,6, 12. All formulae for the one-loop 



4^' {xi,X2) =nfN{ -2 ^^ ^ ^ Ri(xi,X2) + -Y 



Inixif 



+ 



(xi +X2 

1 (-6+19x2 + 25jci) 

6 (JCI +X2){1 -X\ -X2){\ -X\) 



{\-X\-X2)'^ M^\-X\){\-X\-X2) 

1 [-lx\-^x\X2^x\^ -I-^Ixt) 

in(xi) + — — ttt; — TTTi — : — rT^2 



\ {\Xi{x{) ^\Xi{x2))l,2 1 (2-2.^2+.^! 
+ 4 



4xi(l -X\ -X2){\ -X2 



1 (44x1 -31xi;c2 + 13xi2 + 26;c2-26jc2^) 1 

Is Xi(xi+JC2)(l-JC2)(l-JCi-X2) ^^^^^ ^ n 



12xi(l -xi){l -X2){1 -xi -X2] 
ln(x2)2 

l+xi) 24x1 



+2- 



'13 + 6;ci; 



I-X1-X2) X\ 

24 1 



^0^_3(l+^2) 



R{X\,X2) 



8 
1 



(l-.x:2) (^1+^2) 
12 



+ 3- 



(Xl-|-X2)2 X\{\-X2) X\{\-X\){\-X\-X2 



{X\^X2){\-X\-X2) 12 L X\{X\^X2) 
^^"^""'^ ^[Li2(l-Xi)-Li2(l-X2) 



+ - 



1 1 



1 



+ 3 



1 



lx\V (I-JC2) {\-X\){\-X\-X2)- 



N V 3x1^X2^ 



Li2(l -X2] 



1 



03(^1 +X2) +Li3(l —xi —X2) + ^ln(xi +X2)ln(l —x\ — ^2)^ 



1 1 9 o 

--ln(xi)R(jci, 1 -xi -JC2) + -(ln(xi)^-ln(xi)ln(A:2)+ln(x2) )ln(l -xi -X2) 



-^ln(x2)R(x2, 1 -xi -X2] 



+ 



1 



6x\X2 



1^^02(1 -X2) - ^R(^i, 1 -xi -X2) - 2^2 



-f--Li2(l -x\) -Li2(l -x\ -X2) -ln(l -x\ -JC2) -ln(xi +X2)ln(l -xi -X2) 



-^R(X2, 1 -x\ -X2) 



+ 



1 



3X2^ 



1 2 1 2 

- ln(xi ) ln( 1 - XI - X2) - - ln(xi ) ln( 1 - XI - X2) 



1 9 1 

--ln(xi)^ln(l -xi) - -ln(x2)(Li2(l -xi) -Li2(l -X2)) - Li3(xi) -03(1 -xi^ 
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19 

-ln(^i)R(jci, 1 -XI —X2) + — R(xi, 1 -XI -X2) 

o 



+ 



3x1^ 



■ln(x2)ln(l -XI -X2)' 



1 11 
--ln(x2)^ln(l -XI -X2) + -ln(x2)^ln(l -X2) - -ln(xi)(Li2(l - xi) -Li2(l - X2)) 

19 

+Li3(x2) +Li3(l -X2) +ln(x2)R(x2, 1 - x\ -X2) - — R(x2, 1 -x\ -X2) 



^ 1 (19x1 + 13x2) ^^^^ 



+ 



18 XlX2(xi +X2) 

1 (-19 + 19x1+51x2) 
18x2(1 -xi)(l - Xl -X2) 

1 (1-X2 + X1) 

4xi(l -X2)(l -XI -X2) 



x\ -X2) + - 



1 (-1+3X2+Xi) 



In(xi) + 



4X2(1 -^i)(l -x\ -X2) 

1 (1+X1+X2) 



ln(x2)' 



3xi(l -Xl -X2)(xi +X2) 
1 (19-19x2 + 13xi) 



ln(xi)2 

R(X2, 1 -XI -X2) 



1 



+ 



1 



18xi(l -X2)(l -XI -X2) 
2 



ln(x2) + 



lLi2(l-X2) 



Xi(l-X2) X2(l-Xi) Xi(l -Xi)(l -Xl -X2) 



1 

"12 

. l (Li2(l-Xi)-Li2(l 

-Ll2(l -Xl) -Ll2(l -X2)) +- ■ r 

3 Xi(xi+X2) 



3 xi(l -X2) 

(^2 + R(-^l,^2) 

X2) -R(xi,l -Xl -X2)) 



Li2(l -xi; 



(1 -Xl)(l -Xl -X2) 

'3 In(xi) 



1 (Li2(l -Xl) +2R(xi, 1 -Xl -X2)) 
3 



+mn/Ar(- 



1 



X2(l - Xl -X2) 

(2-2x2 + xi) 



(1 -Xl)(l -Xl -X2) 2xi(l -X2)(l -Xl 



2 Li2(l-X2) 

3xi(l -Xl -X2) 

ln(;C2) + ^ R(X1,X2) 



inrif /II 11 

I -— R(xi, 1 -Xl -X2) - 2^^^^2, 1 -^1 -X2) + 



N \2x2^ 

1 (-1+3X2+Xi) 

2x2(1 —x\){\ -Xl -X2) 



X2) 2(l-xi 

1 ln(l — Xl — X2) 

2 X1X2 



-X2) 



In(xi) -- 



1 (I-X2+X1) 



2xi(l -X2)(l -Xl -X2) 



ln(x2) 



(67) 



,(2),fin 
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(xi,X2) =nfN 



1 ln(x2) r^^(l+xi) 59 



36 



50 



36 (1 -X2) L (xi +X2) Xl xi(l -X2) 



+ 12 



(l-3xi)(l-xi) 



1 



+ i^^-^Rl(^l,^2) 



[5 (2-X2) 



18xi(l-X2) (xi+X2)3 



+ ■ 



^1 (2-^2) 



4xi(l-X2)2 

R(xi,X2) 



In(x2)2 + k2 



12xi(l-X2)^ {X1+X2Y 

1 (7xi — 5x2 — 6x1X2 + 5x2^ + xi^) 

12 (1 -X2)xi(l -Xl -X2)(xi +X2) 



Xl(l-Xl-X2) - 

Li2(l -X2) 
In(xi) 



(1 -Xl -X2)(l -X2)(xi +X2) 
5X2^ 



17 



35 



19 1 (1+xi) 1 (l-3xi)(l-xi) 



6^^ 12^^ 12 ' 12(1-X2) 3 (I-X1-X2) 



12x1 



N 
19 



5X2 ^ ^ Xl(l+Xl) 

6x1 (X1+X2) 

1 1 

3xr 
1 



1 (2-X2) 



Li2(l-xi)) 



12xi(l-X2)2 (Xl+X2)2 

"'"7^ (~ ^x~2 ^^^^^ -X2) +ln(l -Xl -X2)^ + ln(xi +X2)ln(l -xi -X2)j 
(X1+X2) 



+ ■ 



18 (1 — X2)xi 3x1 



-ln(xi) ln(l -Xl -X2) + (ln(xi +X2) -ln(x2))ln(l -xi -X2)' 



18 



+ iln(;ci)Li2(l -^i) - ^Li2(l -JC2)ln(jci) -2Li3(l -X2)+2U3{xi +X2) 

+203(1 -XI -X2) -2Li3(jC2) - ^ln(j:2)ln(xi)ln(l -xi -X2) +ln(x2)^ln(l -xi -X2) 

1 5 2 

--ln(xi)R(xi, 1 -XI -X2) - -ln(x2)R(x2, 1 -x\ -X2) -ln(x2) ln(l -X2) 

^ 1 (76x2(1 -X2) +61x1 - 11x1x2) , , . , (17xi^ + 61xiX2 + 38x2^)_ . . . 

^^y^l) H 75 — V ; ^ R(,-*2) i — Xi — X2j 



+ 



36 Xi2(l-X2)2 18xi3(xi+X2) 

1 (-2X2+2X2^-2X1+X1X2) rl , X ,, / n21 , (-«l--«2) 

(-8 + 14x2-6x2^-2xi+xiX2)>, 1 In(xi) l(16xi + 19x2) 



R(xi,X2) +ln(x2) J + 6,^2(,^+,^) R(^^ 1 -^2) 

^ 1 In(xi) , 1(16x1 + 19x2), , 
C2-T^r7^ + o ,.2.,^+,^) ln(l-xi-X2) 

Li2(l -X2) 



12xi2(l-X2)2 12xi(l-X2) 9xi2(xi+X2) 

1 (-10X1^ + 5X1^X2 - 14X1 + 8x1x2 +X1X22 - 10X2 + 10X2^ 



12 (xi +X2)xi2(l —X2Y 

1 (— 2xi'2 +Xi'2x2 + 2X1 — 8x1X2 + 5X1X2^ — 2X2 + 2X2^) 



Li2(l-:ci)) 



12 (xi +X2)xi2(l -X2Y 

J'^^f ( 1 1 /I ^ 1 (-2x2 + 2x2^-2x1 +X1X2) , . x , 

+ ^ b^ln(l -XI -X2) - „ 2n ^2 ln(x2) + 



AT Vxi2 ^ ' 2 Xi2(l-X2)2 ' 2xi(l-X2) 

(xi+X2)^, . A . ^ /I (2-X2) , , , 1 1 



Xi3 ^''^ 



Xl4^^'''"(xi,X2) 



^(^[^ln(l-xi-X2) 



YQg + 4C2 + 2ln(-«i+^2) ln(l-xi-X2) 



1 (52x1^x2 +47x1X2^ - 45x1X2 - 14x2^ + 14x2^ - 19x1^ + 19x1^) 



36 

(xi +2x2) 

6X2 



(1 -Xi -X2)(xi +X2)X2 



R(xi, 1 -xi -X2) 



+ -ln(l -xi)ln(xi)2--ln(x2)Li2(l -X2) + -Li2(l -xi)ln(x2) 



1 2 1 1 

+ -ln(xi)ln(l -xi -X2) +Li3(xi) +Li3(l -xi) + -Li2(l -xi -X2) 

2 Jo 



(X22+Xi2 + 4X1X2) 



12X1X2 

-Li3 ( 1 - XI - X2) - Li3 (xi + X2 



1 1 

+ -ln(x2)ln(xi)ln(l -xi -X2) - -ln(xi +X2)ln(l -xi -X2)^ 



1 (x22 + 5x1 2 + 12x1x2) 



24 



X1X2 



ln(xi)R(xi, 1 -XI -X2) 



1 (-43x1+31x1^ + 63x1X2 + 12-12x2) 
36 (1 — xi)(l — xi — X2) 

1 (2x1^ +X22 + 6x1X2) 



ln(xi) + - 



1 xi(l— xi— 3x2) 



8 (1 -xi)(l -xi -X2) 



In(xi)^ 



4085 



1 



24 X1X2 

1 (xi + 2X1X2 - I+X2) 



In(xi) ln(l-xi-X2)-— — + — Ca 



2592 72 



1 XI (-1+4x2) 



24 (1 -X2)(l -XI -X2) 24 (1 -X2)(l -XI -X2) 



R(xi,X2) 
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+ 



24 



(1-^2) (1 

1 (-82x1 4 



15 

-xi) 



2(4x2-5) 



8X2 



(1-X1-X2) (X1+X2) 

h 14x1X2 + 41 -41x2) 



Li2(l-xi)) 



72 (1 -xi)(l -xi -X2) 

1 (31 —56xi —31x2 + 25x1^ — 20x1x2) , . . 
+^ n ^ ln(xi) + 

36 (1 -xi)(l -XI -X2) 

^(9x1^ + 9x1+58) _ (4xi3- 51x1^ + 58-69x1) 

(1 -XI -X2) 

1 (2-4x1 



ln(xi)^--C3 + 



36 



4345 
2592 



72 



1 



(l--^2) 
(2xi2+XiX2+2x2^^ 



+ 



12 

1 

72 



11 + 



(xi +X2)2 

9 



■Rl(xi,X2) + 



-33x^ + 58x1+58- 

(1 -XI -X2)(xi +X2) 

- 4x2 + 2xi^ + X1X2 + 2x2^' 



58 



(I-X2) 



+ 



45 



12 

9(4x1 + 1) 



(1 -XI -X2)2 

36x1 



(1-xi) 



■R2(xi,X2) 



Li2(l -xi) 



1 xiX2ln(xi) 



(1-Xl) (I-X1-X2) (xi+X2)J 

1 (1340x1 - 1023x1X2 - 670x1^ + 688x1^x2 - 335) 



8(l-xi-X2)2'" 144 (l-^i)(l 

1 (-31 + 13xi)x2 + 31(l-xi)2 



-X2){1 -XI -X2) 



fc) 



+i%nfN 
-21n(l 



24 
72 



(1 -xi)(l -xi -X2) 

1 +xi + 3x2) 



ln(xi) + 



illnf 1 / xi( 



1 X1X2 — 8x1 + 2 + 4x1^ 

8 (1-Xi-X2)2 
1 



R(xi,X2) 



-XI -X2) 



4V(l-xi)(l-xi 

191\ , 
j + (xi ^ 



54 



ln(xi) (x2 + 2xi)R(x2, 1 -xi -X2) 

-xi) Xl 

■X2\ (69) 



(2), fin. ^ ln(xi) 

c\2 (xi,X2) = n/A^( 3 



(xi+X2)2 2xi(l-X2) 



1 T 1 

--ln(x2)^--^2+Li2(l-X2) 

z o 



1 (13x1^ + 36x1-10x1X2-18x2 + 31x2^) 



18 
1 



12x1 (xi +X2Y 
(l+2xi) 



(Xi+X2)2xi(l-X2) 

5x2 + 42x1 + 5- 
6 1 



(1+xi)^ 
(I-X2) 
1 



1 / , , (-X2^-2X1+4X2 + X1^) ^ 

ln(x2) + ^-^^-p Ri(xi,X2) 

_ ^ (1-3x1 + 3x1^) _ xi^ 



+ 



+ 



Xl(xi+X2)^ {Xl^-X2y 12xi(l-X2) 
Hf ( 19 1 



(I-X1-X2) {X\+X2) 
(Li2(l-Xi)-Li2(l-X2))- 



R(xi,X2) 



(xi +X2)xi 



9 xi(xi +X2) 18 



1 (13x1^-29x1+42x1X2-38x2 + 38x2^), , , 
+ — ^ : — ^ ln(x2) 



Xi2(xi +X2)(1 -X2) 



ln(l - Xl -X2) 



l(2-2x2-xi) 2 1 (-32x1-19x2 + 16x1^ + 35x1X2 + 19x2^) 

4 Xi2(l-X2) ^^""'^ ^9 Xi2(xi+X2)2 

(1 - XI -X2) J|in(^2)^- ^ln(xi)ln(x2) + ^ln(xi)2 + ln(xi +X2)ln(l -xi -X2) 



3xr 



■ln(x2)ln(l -Xl -X2)|ln(l -xi -X2) -ln(x2)^ln(l -X2) - ^Li2(l -X2)ln(xi) 
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-2Li3(l -;C2) - 203(^2) +2Li3(l -xi -;C2) +2Li3(jci +X2) + -ln(jci)Li2(l -xi] 



-^ln(xi)R(xi, 1 -x\ -xi) - ^ln(x2)R(^2, 1 -x\ -X2) 



1 

+ - 



In(xi) 



2xi{xi +X2) 



1 {x\^ + 2x\X2+ xj^ — 2xi —X2) 

3 Xl'^{xi+X2)^ 

(2-2x2-xi) 



ln(l —xi —X2)^ + ln{xi +X2)ln(l —xi —X2) 



+Li2(l -xi -X2) 



+ 



12xi2(l-JC2) 



R{xuX2) + 



{4xi +X2+Xi^ —X2^ 



^2 



12x1 



1 



(I-X2) X\ 



6xi2(xi +^2)2 
Li2(l -xi 



Li2(l --^2) 
Ylx\{x\ +JC2) 



Xi(jCi +JC2)2J 12xi(xi+X2) 

6 



R(jci, 1 -^1 -X2) 
6 - 



(I-JC2) JCl (X1+X2) 



+ 



1 R(X2, 1 — -JCl —X2) 
18 Xi(xi+X2) 



(I-JC2) JCi {X\+X2)\ 

17 + 61^-^-38 ^^^^-^^^ . 



, '^"/ / (■^1^ + 2x1-^2 + -^2^-2;ci -X2) , . 1 (2-2;c2--^i) , , . 



1 



X\^{X\ +X2)^ 
( 1 — Xi — ^2 



R{X2: 1 -Xl -X2) 



2 ^l2(l-X2) 

2 ^1 1-X2 



We have introduced the function R{xi,X2), which is well known from ref. 



(70) 



R(-«l,^2) 



^ln(jci)ln(jC2) -ln(xi)ln(l -jci) + ^^2-Li2(xi) ) + {xi ^ X2) . (71) 



In addition, it is convenient, to define the symmetric functions Ri(a:i,a:2) and R2(-^i,-^2) as 



Rl(jCi,JC2) 



R2{xuX2) 



ln(xi)Liij 



Xl \ 1 

— ■ ,xi +X2 - -ln(l -Xl -JC2)C2 -ln(A:i)Li2(;ci +X2) 

Xl +X2 J 2 



4ln(xi)ln(x2)ln(l -xi -X2) -Lii,2 ( ^^^^ — 

2 y^i +^2 



Xl +X2 



-Li2a 



Xl 



Xl +X2' 



Xl +X2 ] +Li3(xi +X2) ] + {Xi ^X2) . 



1 7 3 9 1 

-ln(xi)^ln(l-xi) --ln(xi)^ln(x2) --ln(xi)Li2(l -xi) 
1 



ln(x2)Li2(l -JCi) -Li3(l -xi) +Li3(xi) + {xi ^ X2) 



(72) 



We see in eqs. (]67|)-(|70|) that the multiple poly logarithms all have simple arguments of a partic- 
ular structure, that can easily be continued analytically. Details of this procedure are discussed 
in appendix 
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5 Conclusions 



Determinations of the strong coupling constant from data for e+e 3 jets demand next- 
to-next-to-leading (NNLO) theoretical predictions in perturbative QCD. In this paper we have 
calculated the fermionic contributions to the two-loop amplitude e+e qqg. Furthermore, we 
have obtained the full one-loop amplitude to order £^ in dimensional regularization, needed for 
the interference of the one-loop amplitude with itself. 

We have used a systematic, fast and efficient method for the calculation of loop integrals. Our 
procedure allows for a direct evaluation of all scalar integrals in arbitrary dimensions and with 
arbitrary powers of propagators by means of nested sums. The approach relies on the ability to 
solve all nested sums in terms of multiple poly logarithms, which appear to be the natural class 
of functions for virtual corrections in perturbative QCD. As a consequence, the results allow 
for analytic continuation in a straightforward manner. Therefore, they apply also to (2 + l)-jet 
production in deep-inelastic scattering or to the production of a massive vector boson in hadron- 
hadron collisions. 

The results presented in this paper represent one contribution to the full next-to-next-to-leading 
order calculation of e+e 3 jets. At the same time, it provides an important cross check on the 
results recently obtained by Garland et al. p3|] with a completely independent method. Extend- 



ing our approach to the calculation of the remaining contributions to the two-loop amplitude for 
e+e^ qqg, i.e. the terms which are not enhanced by nf, can be done along the lines of section 

i 
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A One-loop amplitude 

Here we list the result for the one-loop function " to order £?. The remaining functions 
^(i),fin^ . _ 2 4 5^ are of a similar length, and in order to save a few trees, we do not list them 
explicitly here. All results for the one-loop functions cp"*'*^", / = 2,4,6, 12, together with the 

f2) fin 

fermionic contributions c) i = 2,4,6, 12 to the two-loop amplitude can be obtained as a 



FORM file from the preprint server |http://arXiv.or^ by downloading the source file of this article. 



Furthermore they are available upon request from the authors. 



Xl{l—X2) V2 {l—X2)Xl {l—Xl—X2)Xl- 
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XIX2 {l-X2)xi {XI+X2)X2 {xi+X2)^X2 



X\^X2 



^ ln(l-xi-X2)-ln(;c2) ^ r Sy-xx-X2) ^ 

X\^ Xy" 



(X2,l -XI -X2)^ 



(e + nce^ ^ (4-ln(.i)-ln(l-.i-.2)) _ (lzM(in(l -^2) -3)ln(l -x, -X2) 



N 



X\X2 



Xi^X2 



+ 1 (3 -In(x2))ln(x2) - ^i+^^^ R(;ci, 1 -XI -X2) + '^-^\^R{X2, 1 -XI -X2) 
XI ^ X1X2 XI 

1 (61n(x2)-ln(x2)2) _2(2-ln(l -xj -X2)) (4-ln(l -xi -X2)) 



2 (1-X2)xi 
(1 -Xl -X2) 



xr 



(xi + X2 )X2 (xi + X2 ) ^X2 

(31n(l -X2) -ln(l -XI -X2))ln(x2)^-2Ri(l -xi -X2,X2) 

+ (2Li2(l -X2) -4^2-ln(l -XI -X2)^)ln(x2) + 6Li3(l -xi -X2) +6Li3(x2) 
+ ln(xi+X2)ln(l -Xl -X2)^-2(Li2(l -xi -X2) + C2)Ml -xi -X2) ^ 
2^/1 (-2lC2 + 36-91n(x2)+ln(x2)2 



ln(l - Xl -X2) 



■ln(x2) + -: 



6Li3(xi) 



Xl{l-X2) 2xi(l-xi-X2) 
-6Li3(x2) - 31n(x2)2ln(l -X2) + (-31n(l -xi) + ln(x2))ln(xi)2 + 41n(x2)C2 

-(2Li2(l -Xl) -ln(x2)^-4C2)ln(xi) -6R(xi,X2) +2Ri(xi,X2) -2Li2(l -X2)ln(x2) ) 

/ I 
N\6' 



■^1(1 -•^2) ' ■ XiX2 

(8 -7^2 -41n(l -XI -X2) +ln(l -xi -X2f) 1 ln(l -xi -X2) 



(xi +X2)X2 

+61n(l — Xl — X2) — ln(l — Xl — X2)^) + 



(2IC2-24 



2, , (I-X1+X2) 



XlX2^ 



3 (xi +X2)^X2 

Ri(xi, 1 -Xl -X2) -Li2(l -xi)ln(xi) 



1 3 
+ -(ln(xi) -ln(xi +X2))ln(l -xi -X2)^ + 21n(xi)C2 - ^li^C^O^^l -xi) - 3Li3(xi) 

+^(2^2 + 202(1 - Xl -X2) +ln(xi)2)ln(l -XI -X2) - 3Li3(l -xi -X2)] 

Li2(l -X2)ln(x2) -Ri(l -Xl -X2,X2) +3Li3(l -Xl -X2) +3Li3(x2) 



xr 



1 



+^(ln(xi +X2) -ln(x2))ln(l -xi -X2)^-Li2(l -xi -X2)ln(l -xi -X2) 



1 o 

+3R(x2, 1 - Xl -X2) - -ln(x2)lln(l -xi -X2) - 31n(l -X2)) 



+ 



ln(l -XI -X2) 



Xl2x2 



6-7C2 



1 2 3 

+-ln(l -Xl -X2) - -ln(l -Xl -X2) 



+ ■ 



1 



X1X2 



8-7C2-21n(l-xi-X2) + -ln(xi)^ 



+^ln(l -Xl -X2)^ -21n(xi) 



+ ■ 



1 

xi^ 



■ln(l - Xl -X2)^ + 91n(l -xi -X2)C2 



23 



1 , 

--ln(l -XI -X2) +4R(x2, 1 -XI -X2) -61n(l -xi -X2) - 31n(x2)^2 + 61n(x2) 



-^In(x2)2+^ln(x2)3 



+ 



(1 -xi -X2) 



XV 



3 2 

-2Ri(l -XI -X2,X2)ln(l -XI -X2) - -C2 



-(21n(x2)+ln(l -XI -X2) + 2R(xi,X2) +4R(xi, 1 -xi -X2) + 5R(x2, 1 -xi -X2))^2 
+ 12Li4(l -xi) - 14Li4(l -xi -X2) - 14Li4(x2) - 3Li2(l -xi)ln(l -xi -X2)^ 

+(2C2 -ln(xi)ln(l -XI -X2) -Li2(l - xi) + ^ln(l - xi -X2)^-Li2(l -X2))ln(x2)2 
-^ln(xi+X2)ln(l -XI -X2)^ + Li2(l -xi -X2)ln(l -xi -X2)^ + 4Lii3 (1, 1 -xi) 

+ ^(ln(l -XI -X2) -41n(l -X2))ln(x2)2 +Lii,i (^j^, 1 -^1) ln(x2)^ 
+4Li2,2(l, 1 -xi) + 4Li3,i (1, 1 -xi) + 2(Li3(l -xi) +03(1 -xi -X2))ln(l -xi -X2) 
-2(Li2,i (1, 1 -xi) +Lii,2 (1, 1 -^i))(ln(l -XI -X2) +ln(x2)) + [2Li2,i (j^. 1 -^1 
-2Li3(l -xi) + ln(l -XI -X2)ln(xi)^ -21n(l -xi -X2)Li2(l -xi) +2Li3(x2) 
+2Lii,i ^ ^ ^^^^ , 1 - -^1^ ln(l -XI -X2) + 2Lii^2 ^ ^^^^^ ,1-^:^1^ + iln(l -xi -X2)^ 

-31n(xi)ln(l-xi-X2)2 + 2Lii,i,i (^:^^, 1, 1 -xj^] ln(x2) - 6Li3,i (^^i^,l-xi^ 

-2Lii,i,2 



-2Li2,i,i 
-2Lii,i,2 
-6Lii,3 



X2 
1 — Xl 

X2 
1 — Xl 

X2 



,1,1 -Xl 
, 1, 1 -Xl 



1 -Xl -X2 



2Lii 



2,1 



1 — Xl 
X2 



-2Li 



1,2,1 



1 ,1-^1 

1 — Xl 
1 -Xl -X2 
1 — Xl 



■2Li2,i,i 
, 1, 1 -Xl 



1 — Xl 
1 -Xl — X2 
1 — Xl 



,1,1-xiJ -6Li2,2 
, 1, 1 -Xl ) -6Li2,2 



, 1, 1 -Xl -6Lii 3 



6Li3 1 



1 - Xl - X2 
1 — Xl 



1-^1 + 



X2 
1 — Xl 
1 -Xl -X2 
1 — Xl 
1 -Xl — X2 
1 — Xl 

ln(l -XI -X2)^ 



, 1 -Xl 

, 1 -Xl 
, 1 -Xl 



+ ln(xi)^ -2Li2(l - Xl -X2) -21n(xi +X2)ln(l -xi -X2) + 6Li2(l -xi) ^2 



+3Lii 



1 — Vi — Xl \ T 

— ^,1-^1 ln(l-xi-X2)^ + 2 

1 -Xl J 



Lii,i,i 



1 -A i -X2 
1 — Xl 



1,1 -Xl 



-Li2,i (^Y^^' ^ ~^^) (tZ^' ^ ^^^^ ) 



(73) 
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B Analytic continuation 



The multiple polylogarithms Li,„^^^ ..^mi (3^^, are analytic functions in k complex variables 

jj, j = l,...,k. To discuss the branch cuts it is convenient to change the variables according to 

Zj = l-yiy2...yj: (74) 

and one has 

Lim^,...,OTi(j,t,---,Ji) = Li„ ( — ^,1-ziJ. (75) 

The multiple polylogarithms have a representation as nested sums. From this represenation one 
deduces that the multiple polylogarithms are real, if all Zj are in the intervall < Zj < 2 and 
(mi,zi) 7^ (1,0). If zi, Zj-i, Zj+i, Zk are fixed in this interval, a given multiple polylog- 
arithm is an analytic function in the remaining variable zj with a branch cut along the negative 
real axis starting from 0. We denote by Re^^. and Im^^. the real and imaginary part with respect to 
the variable zj. In the calculation presented here, the z/s are ratios of two kinematical invariants: 

Zj = ^ (76) 

In fact, the only ratios which occur are: 

-Sl2 1 ~S23-Sl3 

Xi = , l-Xi = , 

—■^123 —■^123 

-*23 1 -^12-*13 

X2 = , l-X2 = , 

—■^123 — ■5'123 

1 "■^IS , -*12-^23 

l-Xi-X2 = , Xi+X2 = . (77) 

-^123 -■^123 

The real and imaginary part of the logarithm is given by 
Rej;Lii(l-z) = ln|z|, 

Im,Lii(l-z) = -imln(^^^^^ =7r[0(5)-0(O]. (78) 

Eq. ( [78| ) defines the sign of the imaginary part for a ratio of two invariants. All imaginary parts 
of higher multiple polylogarithms can be related to the imaginary part of the logarithm. They are 
easily obtained from the iterated integral represantation 

^hni;,...,mi{-^k'! ■ ■ ■ T-^l^ 

" '* * 'f/dt Y'^' <tl J/dt dt 
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and the fact that 



1 



l-z±/0 

The imaginary part in the variable Zj is given by 



Llm;t,...,mi , ^ 7 - --7 , 



1-zi 

mi — 1 



l-z 



T i^r 



dt 



dt 



Q{t-l) 



o ... o 



dt 



ntk-l 



(80) 



(81) 



Using partial integration, the iterated integral is then related to multiple polylogarithms of re- 
duced weight. For the dilogarithm and the trilogarithm one obtains the well-known formulae 



ImLi2(l-z) 
ImLi3(l-z) 



ln(l-z)ImLii(l-z), 
^ln2(l-z)ImLii(l-z). 



(82) 



In the two-loop amplitude there are also multiple polylogarithms depending on two variables in 
the formLia^((l — Z2)/(l — zi), 1 —zi), where (a^b) — (1, 1), (1,2) or (2, 1). In general we have 

Re^i Rez2 Liaz, ( ^ ~ ^i) + ^' ^^^i ^"^^2 ^^ab ( JT^' ^ ~ ^i) 
+i hn^j Re^2 Liafe (7 — —, 1 - Zi J - Im^i Im^^ Liafc (] — —, 1 - Zi 



l-zi 



l-zi 



For the imaginary parts we find: 

Im^i Lin ( ] — 1 -zi 



Im^, Lii2 



l-Zi 

1-zi' 



Lil ( \ — — ) Im^j Lii (1 -z\) , 



l-zi 



1 -zi = -Lii(zi)Lii 



1-Z2 

l-zi 



Im^i Lii(l-zi), 



Im^iLi2i ( ^ — —,i-zi 
* 1 - zi 



Li2 ( ) Im^, Lil (1 -Zi) ■ 



(83) 



(84) 



Im^jLiii ( ^ — —,\-z\ 
1 -zi 



Lil (1 -zi) -Lil 



1-zi 

1-Z2 



Imz2Lii(l-z2), 
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Imz2 Lii2 



Li2( 



Im„ Li2i 



-Z2 



zi)-Li2[|-|^+Lii(z2)Lii Q 



•Z2 



1-zi 



Li2(l -z\) -Li2 



1-Z2 









-zij _ 



Imz2Lii(l-Z2), 



+ Lii(z2)Lii(l-zi) 



Im,2Lii(l-Z2). (85) 



Im^j IniQ Liu 
Imj;j Im22 Lii2 
Im^i Im^j Li2i 



-Z2 



■Z2 



Zl 
■Z2 



■z\ 



,1 -zij = e(zi -Z2) Inizi Lii (1 -Zl) Im22 Lii (1 -Z2), 
, 1 - Zi^ = -Z2)Lii(zi) Im^j Lii (1 - Zl) Imz2 Lii (1 - ^2) 

1-zi^ 



e(zi -Z2) [Lii(zi) -Lii(z2)] Im^-j Lii (1 -zi) Im^z Lii (1 - 22). 



(86) 



The imaginary parts of the harmonic polylogarithms Lia^(l, 1 — zi) can be obtained from the 
formulae above by setting Z2 = zi. Special care has to be taken for the double imaginary part. 
Here one uses 



e(x-l)^e(x-l) 
ox 

and a factor 1/2 appears in the final formulae: 



lim Im^j Im^2 Lin (7 — 1 -z\ 

zi^zi V 1 — Zl 

lim Im^j Imj.2 Lii2 (\ — 1 -Zi 

zi^zi \i-—Z\ 

lim Im^j Im^j Li2i (- — 1 -Zi 

Z2-^Z\ \\—Z\ 



i|,[e(.-i)f 



^[Im,iLii(l-zi)]2, 



1, 



(87) 



--Lil(zl)[Im,^Lil(l-Zl)]^ 



0. 



(88) 



In the one-loop amplitude we encounter additional multiple polylogarithms. We discuss here as 
an example the imaginary parts of Lim ((1 Z3)/(1 Z2), (1 -Z2)/(l - zi), 1 -zi): 

ImziLiiii {\ — — ^,1-Zi^ =Liii flT^'YT^) I^zi Lii (1 -Zl) , 

1-Z3 



I-Z2 l-zi 
Im^2 Lim ( \ — \ — 1 -zi ) = Lii 

' 1 — Z2 1 — Zl 

X Im22Lii(l-Z2), 

^ ^ . , 1 -Z3 1 -Z2 1 
Im^jLiiii ( ,1-Zi 



1-Z2 



Lii(l-zi)-Lii 



(1 






-Z2/. 



1 - Z2' 1 - Zl 



Liu ( Y3^' ^ 1 
11 



1 - Z2 1 - Zl 
1 - Zl' 1 - Z3 



Im^j Im22Liiii 



-Lil 

1 - Z3 1-Z2 



Lil (1 -zi) -Lil 



1-zi 

1-Z3 



Im^, Im^, Lil 11 



lm^2 Ii^zT Lil 11 



1-Z2 1-21 

X lm^2 Lil (1-22) 

1 - Z3 1 - Z2 



1-zi =e(zi-Z2)Lii 



1 ,1 ,1-zi =9(21-23) 

1-Z2 1-21 / 

X Imj.j Lii(l-zi) Imj.3 Lil (1 -Z3) 

1 - Z3 1 - Z2 



Li, 



1-23 
1-Z2 

1-Z2 



Imz3 Lil (1 -23), 
Im^j Lil (1 -zi) 



Li, 



1-22 
1-23 



Im^, Im^, Im,, Lim 



1 , , , 1 - 21 = 0(z2-23) 

1-Z2 1-21 / 

X Im^j Lil (1 -22) Imz3 Lii (1 - Z3) 
I-Z3 1-Z2 



1 — Z2' 1—21 



1-21 

Lil (1 -zi) -Lil 
1-21 ) =0(zi-Z2)e(z2-23)Im,,Lii(l-zi) 





"•)] 




-22/. 



X Imj.2 Lil (1 - 22) Imz3 Lii (1 - Z3) ■ 



(89) 



From these formulae the imaginary parts of Lim ((1 — Z3)/(1 — zi), 1, 1 — zi) can be obtained by 
setting Z2 = 21 . Double imaginary parts in zi and Z2 are given by 



lim Im^, Im^, Lim 



1 - Z3 1 - Z2 



1-21 



lim Im^i Im^, Im^, Lim 



1-Z2 1-21 

1 - Z3 1 - Z2 



1 



rLii 



1-23 
1-21 



[Im,, Lii(l-zi)]^ 



1 - Z2' 1 - 21 



1-21 



1 



= ^^{zi -23) [Im,, Lil (1 -21)]" Im,3 Lil (1 - 23) • 



(90) 



At weight four there are in addition the multiple poly logarithms with indices Li22, Li 13, Lisi, as 
well as Li2ii, Lii2i and Liii2. The imaginary parts of those are obtained in complete analogy. 
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